Local tunneling probe of (110) Yo.gsCao.osE^CuaOr-^ thin films in a magnetic field 
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Scanning tunneling spectroscopy was performed on (llO)-oriented thin films of Ca-overdoped 
Yo.95Cao.o5Ba2Cu307_i at 4.2K, to probe the local evolution of Andreev— Saint- James surface states 
in a c-axis magnetic field. In zero field, we observed conductance spectra with spontaneously-split 
peaks and spectra with unsplit zero-bias peaks. The former showed enhanced splitting with field, 
and the latter showed threshold splitting above finite fields. Although both field evolutions can be 
described in terms of screening and orbital supercurrents, within the framework of d ± ia pairing 
(d=d x 2_ y 2; a=d xy ,s), the enhanced splitting is consistent with only the d+ia state. Our results have 
direct implications on the stability of broken time-reversal symmetry in cuprate superconductors. 



The predominance of d-wave pairing symmetry in 
the high critical-temperature (T c ) superconducting 
cuprates is responsible for a wealth of unconventional 
phenomenal. Of particular interest is the formation 
of Andreev— Saint- James (ASJ) surface states which are 
manifested as zero-bias peaks in the conductance spec- 
tra taken primarily on (110) tunnel junctions 3 *. These 
low-energy ASJ states arise from the constructive inter- 
ference between time-reversed quasiparticles about the d- 
wave line nodes, thus providing direct information about 
the high-T c order parameter (OP)^ 5 -^. In the case of 
YBa2Cu307_5 (YBCO), much attention has been fo- 
cused on the peak splitting observed in zero field at 
low temperatures 3 , primarily on overdoped samples 8 . 
This spontaneous peak splitting has generally been inter- 
preted as evidence for a complex d ± ia OP (d=d x i _ y i; 
a=d xy ,s) with intrinsically broken time-reversal symme- 
try (BTRS) 9 . Alternatively, there have also been theoret- 
ical suggestions that the spontaneous peak splitting can 
arise extrinsically, from either electron-hole asymmetry, 
multiband effects or impurity perturbatio n 10 ' 11 ! 12 . 

The occurrence of intrinsic BTRS in the thermody- 
namic ground state of a cuprate superconductor would 
have profound theoretical implications on the high-T c 
pairing mechanism^ 3 -. For YBCO it is believed that a 
spontaneous orbital supercurrent, associated with an ia 
OP component, should exist on (110) surfaces. However, 
this spontaneous current has not been directly observed. 
In earlier tunneling studies of YBCO, the existence of 
this spontaneous current was inferred through a Doppler 
interpretation of the magnetic-field evolution of the ASJ 
peak splittin g 14 ' 15 . This interpretation assumes that the 
field-induced screening current enhances the peak split- 
ting through Doppler-like energy shifts of the ASJ sur- 
face states. More recent tunneling experiments on YBCO 
have shown that the rate of peak splitting with decreas- 
ing field is invariant as screening is suppressed with film 
thinness, suggesting that a non-Doppler mechanism for 
the field splitting is also at pla y 16 ' 17 ' 18 ' 19 . 

Another puzzling question about the Doppler interpre- 
tion is why, of all the tunneling data reported on YBCO 
to date, only field enhancement of the peak splitting has 



been seen. Since the two possible states of d ± ia cor- 
respond to oppositely orbiting supercurrents, their spec- 
tral manifestations in a c-axis field should be distinctly 
different. Namely, the peak splitting is expected to in- 
crease/decrease if the screening and orbital supercur- 
rents are parallel/antiparallel. Detailed calculations of 
the peak splitting have indicated a variety of field evo- 
lutions, i.e. asymmetries, thresholds and discontinuities, 
depending on the relative strength of a to cJ2°-. Of these 
possibilities, only the threshold effect has been reported, 
and only in underdoped YBCO 2 ^. It should be noted 
that all of these prior in-field measurements were done 
on macroscopic planar junctions, which are not sensitive 
to local variations in the orbital supercurrents. Recent 
Josephson tunneling and magneto-optics experiments on 
Sr2Ru04 have indicated that a complex OP may actually 
be granularized into time-reversed domains, particularly 
in the case of chiral symmetr y 22 ' 23 ' 24 . 

This Letter reports on the first microscopic ex- 
amination of these issues, using scanning tunneling 
spectroscopy (STS) on (llO)-oriented thin films of 
Ca-overdoped Y .95Cao.o5Ba 2 Cu307_5 (Ca-YBCO). By 
virtue of its nanometer junction size, STS provides a di- 
rect probe of the high-T c OP at coherence-length scales. 
In zero field at 4.2K, we observed conductance spectra 
with spontaneously-split peaks as well as spectra with 
unsplit zero-bias peaks. The former spectra showed en- 
hanced splitting in an increasing c-axis field, and the lat- 
ter spectra showed threshold splitting above finite fields. 
We analyze these spectral evolutions versus field as the 
combined effects of screening and orbital supercurrents. 
Although both field evolutions can be described within 
the framework of d ± ia pairing symmetry, the former 
is consistent with only the d+ia state. The apparent 
absence of the d-ia state has direct implications on the 
stability of BTRS in cuprate superconductors. 

The overdoped Ca-YBCO films measured in this work 
were grown by DC sputtering on (llO)-oriented SrTi03 
(STO) substrates, with a buffer layer of PrBa2Cu307_5 
to ensure proper epitaxy. The films were typically 
160nm thick and showed T c rj 85K, consistent with Ca 
overdoping^ 5 .. The films were transported in airtight con- 
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FIG. 1: (color online) (a) Magnetic field evolution of a 
spontaneously-split dl /dV peak spectrum taken by STS with 
a Pt-Ir tip on a (llO)-oriented Yo.95Cao.o5Ba2Cu307_5 thin 
film at 4.2K. The spectra, staggered for clarity, show the peak 
splitting S increasing with the c-axis field H. (b) Schematic of 
the in-field STM orientation, (c) Plot of S versus H. (d) STM 
image of a film showing strongly-oriented surface morphology. 



tainers and reannealed in flowing O2 before being cooled 
with ultra-pure He gas to 4.2K. The STS apparatus was 
specially designed to sit horizontally in a c-axis magnetic 
field, as shown in Figure [T^b) ■ STS measurements were 
made with Pt-Ir tips by suspending the feedback and 
measuring the tunneling current / versus bias voltage V. 
STM setpoints of InA and lOOmV were typically used. 
For each STS junction, fifty I-V sweeps were averaged 
and then differentiated to give the conductance dl/dV 
spectrum at each field. Ramping between field points 
was done very slowly at ~0.19T/min. Since our STM 
was made entirely from non-magnetic materials, the sta- 
bility of our junctions was mainly limited by piezoelectric 
drift, estimated to be less than Inm at 4.2K from STM 
images taken on NbSe2. 

In zero magnetic field, we observed dl /dV spectra with 
spontaneously-split peaks as well as spectra with unsplit 
peaks, varying over distances as short as ^3nm. Figure 
(Ha) shows the field evolution of a spontaneously-split 
peak spectrum. The magnitude of the peak splitting 
6(H), defined relative to zero bias, is seen to increase 
from 1.3mV to 3.3mV as the field H is increased from 
OT to 3.5T. This overall trend of dS/dH > is summa- 
rized in FigQJc). The relative heights of the split peaks 
appear to change with the field; the negative-bias peak 
is taller than the positive-bias peak at low fields, but 
this asymmetry reverses polarity as H reaches « 3T. Al- 



though such peak-height asymmetry was generally seen 
on our films, its bias polarity and field dependence varied 
spatially. Also worth noting in the spectra of FigQJa) 
are the gap-like features denoted by a slope change in 
dl/dV near ±I8mV, which become less distinct from the 
spectral background with increasing field. These gap-like 
structures were also seen in planar-junction studies 2 ^, 
and could be attributed to afr-plane tunneling contri- 
butions associated with surface roughness. Figure Old) 
shows the STM topograph of a typical film surface, with 
similar morphology as seen in previous studies 2 ^. 

Figure [2] shows the field evolution of another 
spontaneously-split peak spectrum, with a larger (5(0) 
and greater height asymmetry. Here the peak splitting 
is also enhanced by the field, as shown by 8(H) plotted 
in the inset, but the polarity of the asymmetry does not 
change with field. It should be noted that 8(0) as large 
as 2mV were seen across our film surfaces. Nevertheless, 
in virtually all junctions showing spontaneous peak split- 
ting, an overall trend of d8/dH>0 was seen, regardless 
of the magnitude of 6(0). For several junctions we also 
noticed small discontinuities in 8(H) at varying levels of 
H. For the spectra shown in FigfJl a discontinuity occurs 
between 1.5T and 2.0T, as is visible in the inset. 

Figure [3] shows the field evolution of a peak spectrum 
with no spontaneous splitting. As field is increased, the 
peak drops in height and broadens in width before split- 
ting to w 0.8mV at a threshold of w 2.5T. This threshold 
behavior is clearly visible in the 8(H) plotted in the in- 
set. The field thresholds we observed showed variation 
across the film surfaces. It is worth noting that the gap- 
like features in these spectra occur near ±10mV in zero 
field and appear to shift asymmetrically in voltage with 
increasing field. Furthermore, the unsplit peaks observed 
in these spectra tend to be taller than the spontaneously- 
split peaks shown in FigQJa) and FigJ^l In general, all 
the spectral peaks we observed are taller than those re- 
ported in planar junctions measurements, consistent with 
the microscopic nature of our STS probe. 

Our STS data are analyzed in the framework of d ± ia 
pairing symmetry with intrinsic BTRS, to elucidate the 
origin of the spontaneous peak splitting. In the pres- 
ence of an ia OP component on the (110) surface, the 
ASJ surface states are expected to shift by ±q about 
zero energy, giving rise to a spontaneous peak splitting 
(5(0)= a 15 ' 27 i 28 . At zero field, the variation in 8(0) seen on 
our films indicates that a varies spatially. Previous tun- 
neling studies have found that 8(0) tends to increase with 
overdoping&2&, thus the spatial variation we observed 
could be due to doping inhomogeneity. Alternatively, 
this variation could be due to electronic inhomogeneity, 
which may arise from proximity to a quantum critical 
point 8,13 . Here it is worth noting that the combination 
of zero-bias peak, smaller gap size and threshold splitting 
seen in Fig [3] resembles the planar-junction data taken on 
underdoped YBCOSI. 

In the present literature on tunneling in YBCO, the 
existence of a spontaneous orbital current associated 
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FIG. 2: (color online) Magnetic field evolution of a 
spontaneously-split dl /dV peak spectrum taken by STS with 
a Pt-Ir tip on a (llO)-oriented Yo.95Ca .o5Ba2Cu307_5 thin 
film at 4.2K. The peak splitting <5 observed at 0T (black) in- 
creases with the c-axis field H: 2T (green), 4T (red) and 6.5T 
(blue). The inset plots 8(H) showing the overall trend. 



with an ia OP component is inferred from the 5(H) 
evolutio n 14 ' 15 . In this Doppler scenario, the screening 
current combines with the orbital supercurrent to shift 
the energies of the ASJ surface states. At low fields, the 
screening current alone is expected to cause a linear shift 
v~F-p'sCOS(j)c~, where vf is the Fermi velocity, p s is the su- 
pcrfluid momentum and <fi c models the tunneling cone 29 . 
A more explicit calculation by Fogelstrom et al. has pre- 
dicted a variety of nonlinear 6(H) behaviors, i.e. field 
asymmetries, thresholds and discontinuities, depending 
on the relative strength of a to These nonlinear be- 
haviors arise as the screening and orbital currents couple 
to modify both the amplitude and phase of ia on the 
(110) surface. The discontinuities and threshold effects 
we observed (Figs. 2 and 3) are consistent with this sce- 
nario. It is worth noting that these nonlinear effects var- 
ied spatially in our measurements, giving evidence for in- 
homogeneities in the OP symmetry. These observations 
attest to the microscopic nature of our STS junctions, 
which can resolve spatial variations that are presumably 
averaged out in the macroscopic planar junctions. 

In principle, the d+ia and d-ia states can be dif- 
ferentiated by considering the relative orientation be- 
tween the screening and orbital supercurrents. As il- 
lustrated in Figure 31 d+ia (d-ia) corresponds to par- 
allel (antiparallel) flow between the field-induced and 
spontaneous supercurrents near the (110) surface. The 
model calculations of Ref.— have shown that the ASJ 
surface states for d+ia and d-ia tend to be Doppler- 
shifted in opposite directions about zero energy, causing 
the spontaneously-split peaks to either open up further 
or close in instead^. Remarkably, all the spontaneously- 
split peaks in our measurements showed enhanced split- 
ting with field (dS/dH>0), consistent with only the 
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FIG. 3: (color online) Magnetic held evolution of an 
originally-unsplit dl/dV peak spectrum taken by STS with 
a Pt-Ir tip on a (llO)-oriented Yo.95Cao.o5Ba2Cu307_i film 
at 4.2K. As the c-axis field H is applied, the zero-bias peak 
dips and broadens before splitting to «0.8mV at a threshold 
~2.5T, as shown by the peak-spitting 5(H) plot in the inset. 

d+ia state. This observation is puzzling, since d+ia 
and d-ia are degenerate in zero field, and thus half of 
spontaneously-split peak spectra should show dS/dH<0. 
In this regard, the calculations of Refi2°- also indicated 
that parallel flow tends to be energetically favorable over 
antiparallel flow, thereby preferring d-ia to "flip" into 
d+ia above some threshold field H* . Thus the field in- 
crements used in our measurements provide an empirical 
upper estimate on H* « 0.5T, which is within the the- 
oretical range of 0.03T— 0.63T estimated for the d + ia 
scenario 20 . In a followup study we will more precisely 
measure H* by applying finer field increments in the 
lower field regime. 

The relative stabilities of the d+ia and d-ia states in 
a magnetic field can also be considered in the context of 
OP domain s 22 ' 23 . In this scenario, the OP is granularized 
into d+ia and d-ia domains, whereby the spontaneous 
orbital supercurrents from opposing domains cancel when 
spatially averaged. As a result, it would be difficult to de- 
tect such spatially-varying orbital supercurrents through 
macroscopic techniques. Because of its small junction 
size, our STS measurement provides a microscopic probe 
of the OP, presumably within single domains. Therefore, 
the H* estimated above can alternatively be viewed as 
an upper limit on the field required to induce reversal 
between d-ia and d+ia domains. 

The spontaneous peak splitting we observed indicates 
that a is w 10% of d on our overdoped Ca-YBCO 
(110) films at zero field. This observation seems ex- 
ceptional since bulk measurements have suggested a to 
be less than w 5% of d throughout the superconducting 
regime in the cuprate phase diagram over temperature 
and dopin g 31 ' 32 ! 33 . However, the possibility of OP do- 
mains was not explicitly considered in the analysis of 
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FIG. 4: (color online) Heuristic illustration of STS on the 
(110) surface of a d ± ia superconductor in a c-axis magnetic 
field. Field-induced screening currents are represented by the 
dashed arrow, and spontaneous orbital currents by the loops, 
(a) For d+ia, screening and orbital currents flow parallel near 
the surface, causing the spontaneously-split peaks to open up 
further with field, (b) For d-ia, screening and orbital currents 
flow antiparallel, causing the peaks to close in with field. 



these bulk experiments, thus a re-examination of this 
possibility would be important. Here we emphasize that 
it is entirely conceivable for d-wave pairing to be robust 
in the bulk and still have an ia OP component near (110) 
surfaces due to their pair-breaking natur o 15 ' 27 ' 28 . 

Our analysis indicates that, in order to attribute the 
spontaneous peak splitting we observed to an intrin- 
sic ia OP component, it would be necessary to in- 
voke field-induced flipping of d-ia to d+ia. Alterna- 
tive theoretical calculations have suggested that spon- 
taneous peak splitting could also arise extrinsically, from 
either electron-hole asymmetry, multiband coupling or 
impurity perturbatio n 10 ! 11 ' 12 . First, electron- hole asym- 
metry, which is inherent in the Fermi-surface topol- 
ogy of cuprates^, can cause retroflected quasielectrons 



and quasiholes to experience different barrier strengths, 
thereby upsetting their constructive interference^ -. Sec- 
ond, there is experimental evidence to suggest that the 
plane and chain bands of YBCO are coupled^, thus 
providing a multiband mechanism for the spontaneous 
peak splitting 11 . Although both of these mechanisms 
are plausible, at present it is unclear whether they can 
also account for the temperature and field dependences 
of the peak splitting 8 -^^. Lastly, numerical calculations 
have shown that impurities can also produce spontaneous 
and field-enhanced peak splitting 12 . However, the 8(H) 
thresholds and discontinuities seen in our measurements 
were not predicted by these calculations, thus the impu- 
rity scenario seems unlikely to explain our data. Further 
theoretical work is needed to elucidate these possibilities. 

In summary, we have performed scanning tunneling 
spectroscopy on (110) Yo.95Cao.osBa2Cu307_5 thin films 
at 4.2K. In zero field, we observed conductance spectra 
with spontaneously-split peaks as well as spectra with un- 
split peaks. The former showed enhanced splitting with 
a c-axis field, and the latter showed threshold splitting 
above finite fields. Although both field evolutions can be 
described in terms of screening and orbital supercurrents, 
within the framework of d ± ia pairing, the enhanced 
splitting is consistent with only the d+ia state. The ap- 
parent absence of the d-ia state suggests that it is either 
energetically unfavorable in a field, or that the sponta- 
neous peak splitting is due to an extrinsic mechanism 
other than intrinsically broken time-reversal symmetry. 

This work was supported by NSERC, CFI/OIT and 
the Canadian Institute for Advanced Research. 
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